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The tunnel current and magnetoresistance (TMR) are investigated 
in magnetic tunnel junctions consisting of a spin-filter tunnel barrier, 
sandwiched between a ferromagnetic (FM) electrode and a nonmag- 
netic (NM) electrode. The investigations are based on the transfer 
matrix method and the free-electron approximation. The numerical 
results show that the spin transport depends on the relative magneti- 
zation orientation of the FM electrode and the spin-filter barrier, such 
that the tunnel current reaches its maximum when the magnetic mo- 
ments of the FM electrode and the magnetic barrier are parallel. It is 
also found that the TMR increases with increasing the applied voltage. 
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I. INTRODUCTION 



Since the observation of tunnel magnetoresistancc (TMR) in magnetic tunnel 
junctions [1], there has been a lot of interest in structures consisting of two ferromag- 
netic (FM) electrodes separated by one insulator or semiconductor (FM/I(S)/FM). 
This is due to the possible applications as magnetic field sensors and memory cells 
in magnetic random access memories [2,3]. The TMR in the magnetic tunnel junc- 
tions is limited by the spin polarizations of the FM electrodes. Recently, large 
magnetoresistance ratios as high as 30% were reported in Fe/A^Os/Fe [4] and 



The TMR of such FM/I(S)/FM tunnel junctions can be understood in terms of a 
two-band model in which the energy bands of the FM electrodes are split into spin- 
up and spin-down bands with different density of states at the Fermi energy. When 
the magnetization of the electrodes is parallel, the spin-up (spin-down) electrons 
tunnel from a majority (minority) spin to a majority (minority) spin band, whereas 
in the antiparallcl alignment, the spin- up (spin-down) electrons are forced to tunnel 
from a majority (minority) spin to a minority (majority) spin band. This gives rise 
to a change in the conductance when the magnetizations are switched. During the 
last 10 years, a lot of theoretical papers were published on this topic (see reviews 
in Refs. [6,7]). 

Using magnetic insulating layers, instead of having FM electrodes, one can obtain 
very high values for the TMR [8-10]. Tunneling experiments using ferromagnetic 
semiconductors (FMSs), in particular the Eu chalcogenidcs, as the magnetic tunnel 
barriers [11,12], have shown spin polarization in the tunnel current which in favor- 
able cases exceeds 99% [13]. When a magnetic barrier, which acts as spin filter, is 
used in a tunnel junction, due to the spin splitting of its conduction band below Tc 
(Curie temperature), tunneling electrons see a spin dependent barrier height. In 
this case the probability of tunneling for one spin channel will be much larger than 
the other, and a highly spin polarized current may result. More recently, LeClair et 
al. [9] using a spin-filter barrier and an FM electrode, obtained a large TMR in an 
Al/EuS/Gd tunnel junction which is a new method for injecting spins into semicon- 
ductors. The spin filtering and the magnetoresistance have also been investigated 
theoretically in semimagnetic semiconductors [14]. 

In this paper, using the transfer matrix method, we study theoretically the effect 
of the thickness of the magnetic barrier and the applied bias on the TMR for 
tunneling through a FM/FMS/NM tunnel junction with an asymmetrical barrier. 
We assume that the electron wave vector parallel to the interfaces and the spin 
direction of the electron are conserved in the tunneling process through the whole 
system. 

In section 2, we describe the model and present a general formula for the tunneling 
current through the magnetic tunnel junction. Numerical results for a typical tunnel 
junction are presented in section 3. A brief summary is given in section 4. 



Consider a FM/FMS/NM sandwiched structure in the presence of DC bias V a as 
shown in Fig. 1, where </>l and <pn are the barrier heights in the left- and right-hand 
side of the FMS layer above Tc, respectively. For a tunnel junction with different 
electrode materials, the difference in barrier height at both metal-insulator interfaces 
makes the barrier oblique, resulting in an asymmetric current-voltage behavior. In 
this case, an asymmetry parameter A</> = <pL — 4>r is introduced at zero bias to 
account for the tilted barrier potential [15]. Here we consider the case that the 




II. MODEL AND FORMALISM 
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tunneling electron with energy E x is incident from the left and transmits to the 
right along the x direction. In a frcc-clcctron approximation of the spin-polarized 
conduction electrons, the longitudinal part of the effective one-electron Hamiltonian 
may be written as 

Hx = -^^ + Uj{x) + Va - HX) • a ' (1) 

where m* (j=l-3) are the electron effective masses in three regions labeled in Fig. 
1, and 

( 0, x < , 

Uj{x) = 1 E FL + c/> L - (A<f> + eV a )x/d, < x < d , (2) 
[ -eV a , x > d , 

where E F l is the Fermi energy of the FM electrode and d is the barrier width. V a 
which is a spin-dependent potential, denotes the s—f exchange coupling between the 
spin of tunneling electrons and the localized / spins in the magnetic barrier. This 
term, within the mean field approximation, is proportional to the thermal average 
of the / spins, (S z ) (a 7/2 Brillouin function), and can be written as V a = —I<j(S z ). 
Here, a = ±1 which corresponds to a =|, J., respectively and / is the s — f exchange 
constant in the FMS layer. —h(x) ■ a is the internal exchange energy where h(x) 
is the molecular field in the FM electrode and |h| = ho. Although the transverse 
momentum 7ik|| is omitted from the above notations, the summation over k|| is 
carried out in our calculations. 

The Schrodinger equation for a biased barrier layer can be simplified by a coor- 
dinate transformation whose solution is the linear combination of the Airy function 
Ai[p(x)] and its complement Bi[p(a;)] [16]. Considering all three regions of the 
FM/FMS/NM junction shown in Fig. 1, the eigenfunctions of the Hamiltonian (1) 
with eigenvalue E x have the following forms 

A la e ikl " x + B la e- ikl ' x , x<0, 
ip jtr {x) = { A 2a Ai[ P<J (x)} + B 2a Bi[ Pa (x)}, < x < d , (3) 
A 3a e ik3X + B 3a e~ ik3X , x>d, 



where 



and 



kia = ^2ml(E x + h a)/h , (4) 



2m* 3 (E x + eV a + E FR -E FL )/h , (5) 

are the electron wave vectors along the x axis. Aj a and Bj a are constants to be 
determined from the boundary conditions, while 

X 

P<y(x) = t- , (6) 



with 

A = - 



A 



h 2 d 



2m*(A<j) + eV a ) 



1/3 

(7) 



_ [E x E FL -d> L - V°}d 
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Upon applying the boundary conditions such that the wave functions and their 
first derivatives are matched at each interface point Xj, i.e., ipj,a{xj) = ipj+i^{xj) 
and (m*)~ 1 [d'ipj t<7 (xj)/dx] — (m* +1 )~ 1 [dipj+i,o-{xj) /dx], weobtain a matrix formula 
that connects the coefficients A\ a and B\ a with the coefficients A 3a and B 3cr as 
follows: 



M, 



total 



A 3 a 



(9) 



where 



Mi 



total 



Ai[ Pa (x = 0)} Bi[/v(x = 0)] 
Ai'[ Pa (x = 0)} Bi'[ P<7 (x = 0)} 



Ai[p a (x = d)] Bi[ Pa (x = d)} 
Ai'{p a (x = d)} Bi'[ PlJ {x = d)\ 



-ik%d 

o 



o 



! -1 



ik 3 
ik 3 



1 m 3 

Xq 

1 m t 
A m* 



(10) 



Since there is no reflection in region 3, the coefficient B 3a in Eq. (3) is zero and the 
transmission coefficient of the spin a electron which is defined as the ratio of the 
transmitted flux to the incident flux can be written as 



T a (E x ,V a ) 



k la ml 



1 



M 11 

m total 



(11) 



where M^ tal is the left-upper element of the matrix M to tai defined in Eq. (10). 

At T = K, the spin-dependent current density for the magnetic tunnel junction 
in the free-electron model is given by the formula [17] 



J* 



47T 2 r 



eV a 



E F -eV a 



T a (E x ,V a )dE x 



(E F - E x )T a {E x ,V a )dE x 



EF-eV a 



(12) 



where Eq is the lowest possible energy that will allow transmission and is given 

by Eq =max{-/i , —{eV a + Efr — Efl)} for spin-up and Eq = ho for spin-down 
electrons. 

The tunnel conductance per unit area is given by G — JajVa- In this case, 
the TMR can be described quantitatively by the relative conductance change as 



TMR = 



Gff — G 



G 



TT 



(13) 



where Gjj and G-\[ correspond to the conductances in the parallel and antiparallcl 
alignments of the magnetizations, respectively. 



III. NUMERICAL RESULTS 



Taking the Fe/EuS/Al tunnel junction as an example, we calculate the tunnel 
current and TMR of the junction according to the Eqs. (11)-(13). We have chosen 
Fe and EuS because they have cubic structures and the lattice mismatch is only 
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4% [18]. The appropriate parameters for EuS which have been used in this article 
are: 5=7/2, 7=0.1 eV [19]. The parameters Efl and ho for Fe layer are taken 
corresponding to kp^=1.09 A -1 and fcp|=0.42 A -1 (for itinerant d electrons) [20]. 
In the Al layer Efr — 11.7 eV [21]. The barrier heights at the interfaces, which can 
be derived from the work functions [22] of Fe and Al and the electron affinity of EuS 
[11], are taken as 0^=1.94 eV and </>#=1.7 eV. In practice, the effective masses of 
electrons may differ from that of free electron, but here for simplicity, we assume all 
electrons have the same mass, m as free electrons. We show the numerical results at 
T=0 K. Thus, in the case of parallel (antiparallcl) alignment (S z ) — S ((S z ) = —S). 

In Fig. 2, the TMR is shown as a function of the width of the magnetic barrier 
d, when the bias voltage 14=0.5 V is applied to the junction. At zero temperature 
and for the nearly normal incidence, the electrons with E x near Ep carry most 
of the current. Furthermore, the Fermi energy and the barrier heights at fixed 
temperature, are constant. Therefore, with increasing the thickness of the EuS 
layer from zero, when the thickness approaches half the wavelength of the electron 
wave in the magnetic barrier, the tunnel conductance with the parallel alignment 
increases faster than with the antiparallcl alignment. In this case, wc can expect a 
peak in the TMR. By increasing the thickness further, the tunnel conductance and 
hence the TMR decreases exponentially. 

Figure 3 shows the spin-dependent current densities for spin-up and spin-down 
electrons as a function of applied bias. In the parallel alignment, the tunnel current 
for spin-up electrons is much higher than the spin-down ones, while in the antiparal- 
lel alignment, the discrepancy between the tunnel currents of the two spin channels 
increases very slowly. The origin of this effect is that the tunnel current for each 
spin channel depends on the density of states in the Fe layer and the heights of the 
tunnel barrier. In the parallel (antiparallcl) alignment, the majority spin electrons 
tunnel through a barrier with low (high) height, whereas the minority spin electrons 
tunnel through a barrier with high (low) height, thus we can expect a high (low) 
spin-polarized current. 

It is clearly observed that for both the parallel and antiparallel alignments at low 
voltages the tunnel current curves vary linearly. With increasing the bias voltage, 
the effective width of the barrier becomes narrower and hence the slope of the curves 
increases. This effect is very strong for spin-up electrons in the parallel alignment, 
and causes the TMR increases with increasing the bias voltage. 

In order to reveal another aspect of the spin filtering phenomenon in the system, 
we have displayed the voltage dependence of the TMR in Fig. 4. As the figure shows 
for d=7 and 10 A, with increasing the bias voltage the TMR increases linearly and 
decreases at high voltages. This decrease which is not identical for different barrier 
thicknesses is a consequence of the narrowing of the barrier width. At low voltages, 
the tunneling probability for one spin channel is higher than the other, whereas 
at high voltages, by narrowing the barrier, this tunneling probability increases for 
both channels; thus, the tunnel conductance ratio, Gfi/Gft, tends toward unity, 
and hence the TMR decreases. This dependence of the TMR on bias voltage is 
not similar to the results of the FM/I(S)/FM tunnel junctions [23], because in the 
present study, we have a spin-dependent barrier. In FM/I(S)/FM junctions, both 
the majority and minority electrons see identical barrier heights. In this case, as the 
bias voltage increases, the imbalance between the number of majority and minority 
tunneling states diminishes and, therefore, the TMR decreases. 

From these results one can see that, for the FM/FMS/NM tunnel junctions, the 
spin transport and hence the TMR can be controlled by the applied bias and the 
thickness of the FMS layer. 
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IV. CONCLUDING REMARKS 



In summary, we have shown how using a magnetic barrier and a FM electrode one 
leads to a new magnetoresistive device. The obtained results for spin-filter mag- 
netoresistance show that the spin currents depends on the relative magnetization 
orientation of the FM electrode and the spin-filter barrier. It was also found that 
the low spin-polarized current will be amplified through the magnetic barrier; thus, 
a very large spin-polarized current from the FM electrode is injected into the NM 
electrode and a high value for the TMR is then obtained. 

Our analysis of the TMR can be potentially useful to achieve larger TMR by 
optimally adjusting the material parameters. The two materials systems, ferrites 
and garnets [24] are both insulating and ferrimagnetic above room temperature. 
Therefore, they are well suited for obtaining large TMR in FM/FMS/NM. How- 
ever, for successful spintronics applications, future efforts will have to concentrate 
on fabricating FM semiconductors in which ferromagnetism will persist at higher 
temperatures. 
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Figure 1 A. Saffarzadeh 

FIG. 1. Spin-dependent potential profile for FM/FMS/NM tunnel junctions in the pres- 
ence of a positive bias V a - In the FMS layer, the dashed line represents the bottom of the 
conduction band at T > T c and the thin arrows indicate the bottom of the conduction 
band for spin-up and spin-down electrons at T < Tc- The direction of magnetization in 
the FM layer is fixed in the +z, while the magnetization in the magnetic barrier is free to 
be flipped into either the +z or — z direction, as indicated by thick arrows. The zero of 
energy is taken at the middle of bottoms for majority-spin band and minority-spin one in 
the FM layer. 
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Figure 2. A Saffarzadeh 

FIG. 2. Dependence of the TMR on the thickness of EuS layer. 
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Figure 3. A. Saffarzadeh 

FIG. 3. Dependence of the current densities J CT on bias voltage calculated for d~7 A in 
the parallel (P) and antiparallel alignments (AP). 



10 



50 




20 L- i l i l i l i 

0.0 0.5 1.0 1.5 2.0 

V(volts) 



Figure 4. A. Saffarzadeh 

FIG. 4. Dependence of the TMR on bias voltage calculated for d=7 and 10 A. 
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